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Abstract 

Thermodynamic theory is developed for the ferroelectric phase transition of an ultrathin film in equi- 
librium with a chemical environment that supplies ionic species to compensate its surface. Equations of 
state and free energy expressions are developed based on Landau-Ginzburg-Devonshire theory, using elec- 
trochemical equilibria to provide ionic compensation boundary conditions. Calculations are presented for 
a monodomain PbTiOs (001) film coherently strained to SrTiOs with its exposed surface and its electron- 
ically conducting bottom electrode in equilibrium with a controlled oxygen partial pressure. The stability 
and metastability boundaries of phases of different polarization are determined as a function of temperature, 
oxygen partial pressure, and film thickness. Phase diagrams showing polarization and internal electric field 
are presented. At temperatures below a thickness-dependent Curie point, high or low oxygen partial pres- 
sure stabilizes positive or negative polarization, respectively. Results are compared to the standard cases 
of electronic compensation controlled by either an applied voltage or charge across two electrodes. Ionic 
surface compensation through chemical equilibrium with an environment introduces new features into the 
phase diagram. In ultrathin films, a stable nonpolar phase can occur between the positive and negative polar 
phases when varying the external chemical potential at fixed temperature, under conditions where charged 
surface species are not present in sufficient concentration to stabihze a polar phase. 

PACS numbers: 77.80.bn, 64.70.Nd, 68.43.-h, 77.84.Cg 
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I. INTRODUCTION 



The equilibrium polarization structure of an ultrathin ferroelectric film is strongly affected by 
the nature of the charge compensation of its interfaces. When there is insufficient free charge at the 
interfaces, a locally polar state can be stabilized by formation of equilibrium 180° stripe domains 
that reduce the depolarizing field energy.'^ When electrodes are present, electronic charge at the 
interfaces can stabilize a monodomain polar state, provided that the effective screening length in 
the electrodes is sufficiently small compared with the film thickness. In both cases, the Curie 
point Tc is expected to be increasingly suppressed as film thickness decreases because of the 
residual depolarizing field energy. 

Even when the surface electrode is missing, it has been experimentally observed that a mon- 
odomain polar state can be stable in ultrathin ferroelectric films .'^^'^ This has been attributed to the 
presence of ionic species at the surface that provide charge compensation and reduce the depolariz- 
ing field energy,'^ similar to the adsorbates observed on bulk ferroelectric surfaces.'^ Furthermore, 
experiments have shown that the sign of the polarization can be inverted by changing the chemical 
environment in equilibrium with the surfaccl^^EL Recentl)^^ it was found that when the polariza- 
tion is inverted by changing the external chemical potential, switching can occur without domain 
formation and at an internal field reaching the intrinsic coercive field for certain ranges of film 
thickness and temperature. Thus, through either kinetic suppression of domain nucleation, or the 
structure of the equilibrium phase diagram, an instability point of the initial polar state can be 
reached. This is in sharp contrast to switching by applied field across electrodes, where the con- 
sensus has been that polarization inversion occurs only by domain nucleation and growth at fields 
well below the instability.!^ 

These studies motivate the need to understand the polarization phase diagrams and metastability 
limits for ultrathin ferroelectric films with ionic surface compensation, in chemical equilibrium 
with their environment. While the energy and structure of ferroelectric surfaces compensated by 
ions have been predicted by ab initio c alculations J ^ ^l^"!^"^ !!!! these zero-temperature results have 
been extrapolated to experimental temperatures using simple entropy estimates, and to date have 
not included the effects of interaction with the ferroelectric phase transition and Tc- Here we 
develop a thermodynamic theory for the ferroelectric phase transition of an ultrathin film in an 
environment that supplies ionic species to compensate the polarization discontinuity at the surface 
of the ferroelectric. We use an approach based on Landau-Ginzburg-Devonshire (LGD) theory 
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for the ferroelectric material, with boundary conditions that include both the depolarizing field 
effects that arise in ultrathin films and the creation of ionic surface charge through electrochemical 
equilibria. This new chemical boundary condition is based on a Langmuir adsorption isotherm 
for ions We develop an expression for the free energy of the system and use it to determine the 
equilibrium monodomain polarization states and their stability. For simplicity we do not include 
additional "intrinsic" surface effects or polarization gradients in the ferroelectric.'^ We compare 
and contrast our model for ionic compensating charge controlled by an applied chemical potential 
with existing models for electronic compensating charge controlled by either an applied voltage 
or fixed charge, to elucidate how the present predictions for ionic compensation differ from prior 
work. 

We find that the equilibrium phase diagram of a monodomain ferroelectric film as a function of 
temperature and chemical potential can have a different form than the standard phase diagrams as 
a function of temperature and applied voltage or charge. We present calculations for PbTiOs (001) 
films with a conducting bottom electrode (e.g., SrRuOs), coherently strained to SrTiOa, and with a 
surface compensated by excess or missing oxygen ions.E^^^Slpor sufficiently thin films, we find that 
a nonpolar state becomes stable between the positive and negative polar states, within the range 
of external oxygen partial pressures where there is insufficient surface charge to stabilize a polar 
state. Under these conditions the Curie temperature depends strongly on the oxygen chemical 
potential. 

II. THERMODYNAMIC MODEL 

In this section we establish the electrostatic boundary condition, the ferroelectric constitutive 
relation, and the free energy expressions used to describe an ultrathin ferroelectric film. We con- 
sider a uniformly polarized, monodomain film with uniaxial polarization oriented out-of-plane 
(normal to the interfaces). This should apply to systems such as PbTiOa (001) coherently strained 
to SrTiOs, since LGD theorjl^^ predicts that compressive in-plane strain favors this "c domain" 
polarization orientation. Even if out-of-plane polarization is suppressed by depolarization field 
effects, in this system the in-plane "a domain" polarization orientation is less stable than the non- 
polar phaseP^at temperatures above 360 K. For this case all fields can be specified by scalars since 
their in-plane components are zero. 

To include the effects of incompletely neutralized depolarizing field, we use the simple electro- 
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static model illustrated in Fig. [T| The spatial separation between the compensating free charge in 
the electrodes and the bound charge at the outer surfaces of the ferroelectric leads to residual depo- 
larizing field in the film even when the external voltage Vex is zero (i.e. short-circuit conditions). 
Figure [T] shows the polarization P and displacement D ma ferroelectric film of thickness t sand- 
wiched on the top and bottom by planes of compensating free charge of density ±cr, at a distance 
A > outside the ferroelectric. The planes of bound and free charge lead to steps in P and D, 
respectively. In Fig. [T| P and D are positive and the free charge on the top electrode a is negative. 
The electric field and potential can be calculated from E = (D — P) / and V0 = —E, where eg 
is the permittivity of free space. The internal field in the ferroelectric film is Ein = — (a + P)/^q, 
while the field just outside the film is E\ = —a/eo. 

In a series of early papers, Batra, Wurfel, and SilvermanP^l showed that the results of a more 
complex model taking into account the space charge distribution and nonzero screening length 
in non-ideal metal electrodes could be reproduced by this simple model in which ideal metal 
electrodes are separated from the ferroelectric by a vacuum gap of width equal to the screening 
length, and all bound and free charges reside at the interfaces. This model has been discussed 
extensiveljJ^'^^'^^'^ and used to parametrize the results of ab initio calculations.'^ An alternative 
description in terms of interfacial capacitanceP^^ is equivalent if the interfacial capacitance per 
unit area is identified with eo/A. Recent calculation ij^^ * '^ ! ^"* have shown that the screening length 
for the electrode material can be generalized to be an effective screening length for a given ferro- 
electric/electrode interface. 

We can relate the external voltage Vex across the structure to a and P by integrating the field to 
give 

eoVex = 2Xa + t{a + P). (1) 
The field in the ferroelectric can then be expressed as a function of P and either a or Vex using 

Eir^ = -(a + P)/eo 

_ Vex + 2AP/60 

2\ + t • 

In the latter expression, the second term in the numerator gives the voltage from the residual 
depolarizing field that is proportional to (and opposes) the film polarization. 

In this simple electrostatic model, we assume that the two interfaces have the same screening 
length and work function. In a polarized material, these quantities can depend upon the polar- 
ization magnitude and orientation with respect to the surface, and differences between the two 
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interfaces may arise even if the electrode materials are identical.'^^'^ Since to first order these ef- 
fects simply add a term to 2A in the numerator of Eq. ([2]), which is already a variable parameter 
in our model, we have neglected these differences. The approximations in this electrostatic model 
are not critical in determining the new behaviors we find below for ionic surface compensation 
(which occur even for A = 0), but do provide simple, analogous electronic compensation models 
for comparison. 

To determine the equilibrium polarization in the film, the values of field and polarization inside 
the ferroelectric must simultaneously satisfy both the electrostatic boundary condition Eq. ([2]) and 
the Ein{P) constitutive relation for the ferroelectric. For PbTiOa this can be written aP 

Eirr = f{P) = 2alP + 4al,P^ + 6«niP', (3) 
where f{P) is the derivative of the bulk LGD free energy density 

f{P)=a;P' + a;,P' + aniP', (4) 
and the coefficients a* are those for a coherently-strained film,El 

~ 2eoC ~ sii + S12' 

«33 = "11 + , > (5) 

-Sll + S12 

where Xm is the epitaxial misfit strain of the zero polarization state, Tq is the temperature at which 
changes sign for Xm = 0, C is the Curie constant, and Qij and Sij are the electrostrictive and 
elastic compliance coefficients, respectively.'^'^^'^ Values of these material parameters for PbTiOs 
are listed in Table |l| The misfit strain Xm has a somewhat temperature-dependent valud^^lgf about 
-0.01 for PbTiOa coherently strained to SrTiOa. While for unstressed bulk PbTiOa the fourth- 
order polarization coefficient an is slightly negative, indicating a weakly first-order transition as 
a function of T at Ein = 0, for coherently strained films the coefficient has a positive value of 
5.0 X 10'' Ym^/C^, indicating that the transition is second order-lSI 

The strain normal to the film can be calculated from the polarization using the expressiorP 

X3 = QnP^ + 2su{xm - Qi2P^)/{sii + S12). (6) 

If the effects of depolarizing field are neglected (i.e. for ideal electrodes with A = 0), the Curie 
temperature is determined by the change in sign of c^g, which gives 

T° = To + ^eoCQuxJ{s,i + s^)- (7) 
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Using the XmiT) appropriate for epitaxially strained PbTiOs on SrTiOs, this gives = 1023 K, 
about 270 K higher than in the Xm = case. 

We can determine which of the equilibrium solutions is stable, metastable, or unstable by con- 
sidering the appropriate free energy of the system. For a closed system (e.g. fixed charge), the 
Helmholtz free energy is minimized at equilibrium. The Helmholtz free energy per unit area A 
can be written aP 



^ = t[fiP) + '4El]+2X'-^El 



t 



f{P) + 



2 



eo 



(8) 



2eo 

where the two terms are for the ferroelectric film and surrounding screening regions. For an open 
system (e.g. fixed potential), the Gibbs free energy is minimized at equilibrium. The Gibbs free 
energy per unit area Q is given b)^' 

g = A- V,,a, (9) 

where the difference between the Gibbs and Helmholtz free energies — Vexcx is the electrical work 
done on the system by the external circuit. This difference is in accord with that in a recent 
derivatioiP^ of the energy functionals to be minimized in first-principles calculations at fixed D 
and fixed E. 



III. FERROELECTRIC FILM WITH ELECTRONIC COMPENSATION 

In this section we present the equations of state and phase diagrams for ferroelectric films 
with electronic compensation under controlled voltage or charge conditions, as background for 
development of theory for ionic compensation. Some of the more subtle differences between fixed 
voltage and fixed charge boundary conditions are highlighted. 

Figure |2] graphically shows the equilibrium polarization and field values that simultaneously 
satisfy the constitutive relation, Eq. ([3]), and either the fixed Vex or the fixed a boundary condi- 
tion. A temperature of 644 K was chosen to match one of the experimental conditions previously 
studied.'^S'^Each line in Fig.|2];a) is the fixed Vex boundary condition from the second equality of 
Eq. ([2]) for a particular value of Vex- The deviation of this line from vertical reflects the nonzero 
value of \/t = 10^'^ used to model the screening length in the electrodes. As Vex is varied, this 
boundary condition translates along the horizontal Ein axis. For \Vex\ less than a certain value. 
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there are three intersections representing equilibrium solutions; at larger \ Vex\, there is only a sin- 
gle solution. The marked intersections correspond to solutions that are not unstable, as described 
below. Each line in Fig. |2];b) is the fixed a boundary condition from the first equality of Eq. (|2]) for 
a particular value of a. These lines are nearly horizontal, showing that the field dependence of P 
at constant a is negligible. As a is varied, this boundary condition translates along the vertical P 
axis. In this fixed, uniform a case, there is a single equilibrium solution at all a and T values. The 
behavior is independent of A and t, and, as described below, the equilibrium solution is always 
stable. 



A. Phase diagram for controlled Vex 

At constant V^x, the equilibrium polarization value is that which minimizes Q. Using Eq. ([T]) to 
eliminate a gives an expression for Q in terms of V^x and P, 

G = tf(P)- + + ^^^^ (10) 

Figure |3] shows ^ as a function of P and Vex corresponding to Fig.|2];a). The equilibrium polar- 
ization Peg can be determined by setting the first derivative of Q at constant Vex to zero, giving the 
equation of state 



t dP 



This agrees with the simultaneous solution of the constitutive relation and boundary condition 
shown above, Eqs. (|2]) and Q. 



The stability of the equilibrium solutions of Eq. ( 1 1 1 is determined by the sign of the second 
derivative of Q at that value of P, 



1 d-'g 



t dP^ 



2 \ 



When this is negative, the solution is unstable; when it is positive, the solution is stable or 
metastable. In particular, when there are three solutions, as shown in Fig. [2]^a), the middle one 
near P = is unstable. The values of P and Vex at the instability (limit of metastability) are 
given by the condition that both the first and second derivatives of Q are zero. At this value of 
P, the P(Ein) curve of the constitutive equation, Eq. ([3]), has the same slope as the constant Vex 
boundary condition line, Eq. (|2]), in Fig. ^a). The value of at the instability is the intrinsic 
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coercive field for the film/electrode system with parameters t and A, taking into account the effect 
of depolarizing field. 



The solution of Eq. ( 12 1 for P = gives the condition for the Curie temperature Tc which can 



differ from the value for A = 0. The change in Tc due to a nonzero screening length is given 

— 2AC 

AT = Tc-T?. = — . (13) 

Because the Curie constant C is much larger than Tq for typical ferroelectrics, stability of the polar 
phase requires A << t. Even a ratio \/t = 0.001 gives AT = —300 K for PbTiOa. Effective 
screening lengths A calculated from first principle^^^^ vary between zero and 0.02 nm for various 
electrode-ferroelectric interfaces . 

Figure |4] shows the equilibrium polarization, strain, and free energy of the stable and metastable 
solutions as a function of Vex- These indicate the possible polarization hysteresis and strain but- 
terfly loops. Two equilibrium solutions (one stable and one metastable) corresponding to oppo- 
sitely polarized states exist when \ Vex\ is smaller than the instability. The stable solution switches 
between positive and negative polarization at V^x = 0. At values of V^x in the metastable re- 
gion between the equilibrium and instability points, polarization switching requires nucleation 
of domains of the opposite polarity. The nucleation barrier becomes zero when Vex reaches the 
instability.'^ At values of Vex exceeding the instability, switching occurs by a continuous process 
without nucleation. 

Figure [5] shows the equilibrium polarization phase diagram as a function of Vex and T. While 
there is a first-order transition phase between positive and negative polar ferroelectric (F-i- and F-) 
phases, this terminates at Tc in a second-order transition to the nonpolar paraelectric (P) phaseP'. 
When T is varied at nonzero values of Vex, there is no phase transition between the nonpolar and 
the stable polar phase. The dashed red and blue curves are the limits of the metastable F-i- and F- 
phases, respectively. The polarity switching transition driven by changing Vex at fixed T requires 
nucleation under conditions inside (below) these curves, and is continuous (non-nucleated) outside 
(above) these curves. 

The nonzero screening length A not only depresses Tc below T^, but also produces inverted 
electric fields in the film in the region near the phase boundary (small \Vex\)- Figure [6] shows 
the internal field as a function of Vex and T in the vicinity of Tc- The inverted electric fields 
extend into the nonpolar phase in the region between Tc and T^. Thus, when a small external 
voltage is applied to a film in this region, the equilibrium field in the film is opposite to the applied 
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field. Close to Tc, the magnitude of this inverted field is larger than that of the applied field, 
producing an (internal) voltage gain in a passive device that diverges as Tc is approached. It 
has been proposed that such "negative capacitance" could be used to improve the performance of 
nanoscale transistors.^'' The conditions for zero field, shown as black dashed curves in Fig.[6[ can 
be obtained from Eq. ^ as Vex = — (2A/eo)Po: where Pq is the zero-field spontaneous polarization 
of the epitaxially strained film."^ Unlike Tc, these boundaries are independent of film thickness. 



B. Phase diagram for controlled a 

An alternative to controlling the voltage Vex across the electrodes is to control the current flow to 
the electrodes, and hence the free charge a. The equilibrium polarization at fixed a is determined 
by minimizing the Helmholtz free energy for the system given by Eq. ([8]). This is plotted as a 
function of P and a in Fig. |7] Setting to zero the derivative of A at constant a gives the equation 
of state 



t dP 



o + P 

f{P) + (14) 
eo 



This agrees with the constitutive relation, Eq. Q, verifying that the correct equilibrium states are 
predicted by this free energy for fixed a. Since C/Tc is very large for typical ferroelectric s, there 
is a only one minimum at each value of a. The equilibrium values of A and Vex are plotted versus 
cr in Fig. [8] To a good approximation, this equilibrium solution corresponds to Peg ~ —a, and 

Aeg ~ tf{a) 



Vel' ^ tf\a) + 2\-. (15) 
eo 

If these approximate expressions were plotted with the exact results in Fig. [8| the curves would be 
indistinguishable. 

The stability of this equilibrium solution with respect to variations in P can be evaluated from 
the sign of the second derivative of A with respect to P, 



1 d^A 



f'iP) + -. (16) 
eo 



t dP^ 

Because C/Tc is very large, the first term is negligible, and the second derivative is always posi- 
tive. Unlike the fixed Vex case, the equilibrium solution Peg ~ —a is never unstable with respect 
to fluctuations in P, and there is no phase transition. This point is often not recognized - if the 
surface charge is controlled and kept uniform, any value of P in the film may be stably formed. 



While at fixed surface charge density a the equilibrium solution is always stable with respect to 
polarization variations, instability can occur with respect to spatial nonuniformity in a. The free 
energy Aeq of Eq. ( 15 1 is a double well, as shown in Fig. [8];b). The minima occur at values amin 



given by solutions of 

= « /V) + ^. (17) 

t da teo 

This gives V^^{amin) = 0. If the controlled parameter is the net charge density on the electrode 
a, then for |a"| < |o"mm| the system can lower its free energy by forming a two-phase mixture of 
domains of opposite polarity. The extent of this two-phase region is shown by the black dotted lines 
in Fig. Isj The local surface charge a = ±\amin\ will have opposite sign for oppositely polarized 



domains. At equilibrium, the fraction of positive domains Xpos will be Xpos = (1 — d-/\amin\) /2. 
The equilibrium value of V^x is zero in this polydomain region and the free energy densities of the 
oppositely polarized domains are equal. 

Here we assume that the in-plane size of the domains is sufficiently large compared with the 
film thickness that we can neglect the excess free energy of the domain walls and the in-plane 
components and nonuniformity of the polarization and field near the domain walls. When there 
is incomplete neutrahzation of the depolarizing field by compensating charge (e.g. when A is not 
zero), the free energy can in some cases be reduced by the formation of equilibrium 180° stripe 
domain^^ with an in-plane size similar to or less than the film thickness. For such fine-scale 
domain structures, the domain wall energies and field variations are not negligible. For simplicity, 
we do not consider these cases here. 

The surface charge density is a conserved order parameter, so instability with respect to spatial 
nonuniformity occurs for magnitudes of a smaller than the spinodal values given by 

0=^^^/» + ^. (18) 
Since A << t, this expression gives the same result as the instability with respect to uniform P 



variations at constant Vex, Eq. (12). 



Figure |9] shows the equilibrium polarization phase diagram as a function of a and T, while 



Fig. 10 shows the internal field in the vicinity of Tc. These exhibit a two-phase field between 
the single-phase positive and negative polar phases (F-i- and F-). The dashed red and blue spin- 
odal curves are the metastability limits of the F-i- and F- phases, respectively. The polarization 
behavior is especially simple, since we have P ~ —a in the single-phase regions, independent 
of T. The suppression of Tc is the same as in the controlled Vex case, and the equihbrium and 
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instability curves correspond exactly. However, transformations driven by controlling either Vex 
or a follow different paths. If V^x is kept constant, the parent phase will remain metastable and 
will be entirely consumed by the stable phase. If a is kept constant, Vgx will decrease to zero as 
the fraction of inverted domains grows, reaching an equilibrium two-phase state. The controlled 
potential and controlled charge phase diagrams for ferroelectrics. Figs. |5] and |9[ are directly anal- 
ogous to controlled chemical potential and controlled composition phase diagrams for an alloy 
or fluid exhibiting phase separation.'^^ In particular, the instability in ferroelectrics is a spinodal 
boundary, and the continuous transition that occurs at fixed Vex in the unstable region is equivalent 
to spinodal decomposition of an alloy held at constant chemical potential. In this case, unlike 
the usual fixed average composition constraint for an alloy, the continuous transition will result in 
a single-phase (monodomain) final state. Spinodal transitions from monodomain to polydomain 
states in ferroelectrics at fixed a have recently been modeled.'^ 

While the conclusion that the P ^ —a solution is always stable with respect to P fluctuations 
for any constant a may seem practically irrelevant for the electronic compensation case where 
the system is unstable with respect to a nonuniformity, in the case of ionic compensation this 
conclusion can be important. As we shall see, for ionic compensation the system can be stable 
against a nonuniformity, and the phase transition to a polar state can be completely suppressed for 
a range of applied chemical potential. 

IV. FERROELECTRIC FILM WITH IONIC SURFACE COMPENSATION 

Now we consider a ferroelectric film without a top electrode, but with its surface exposed 
to a chemical environment that can supply free charge from ionic species. The amount of free 
charge supplied will depend on the chemical composition of the environment and the external 
voltage Vex that it sees on the surface. We will use the same electrostatic boundary condition, 
Eq. ([T]), constitutive relation, Eq. ([3]), and free energy, Eq. ([8]), employed above for the electronic 
compensation cases, treating the ions as residing in a plane at a distance A above the surface. 
Rather than solving for the polarization for a given value of Vex or a, we wish to obtain the 
equilibrium polarization for a given composition of the environment. 

To obtain the relationship between a and Vex due to this chemical equilibrium, we develop an 
expression based on those for adsorption of ions used in electrochemical systems. We treat the 
external chemical environment as an electrolyte that is in contact with both the surface of the film 
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and the bottom electrode (e.g. via pinholes in the film away from the region of interest). In order 
for surface ions from the chemical environment to produce an electric field across the sample, as 
observed in experiments j^^^ the electrons involved in creating the surface ions must have such a 
path to reach the bottom electrode. 

We can write a generalized surface redox reaction between oxygen in the environment and a 
particular surface ion i, 

lonSite + — O2 O -ZjC" + Ion^% (19) 

where Ui is the number of surface ions created per oxygen molecule, and zi is the charge on 
the surface ion. In this formalism, Ui and Zi change sign depending upon whether positively or 
negatively charged surface species are involved. For example, if the surface ion is a doubly- 
negatively-charged single-atom adsorbed oxygen, O^^, so that rij = 2 and Zi = —2, the redox 
reaction is 

Va,+ ^02 + 2e- oO^-, (20) 

while if the surface ion is a doubly -positively-charged single-atom missing surface oxygen, V^"^, 
so that Hi = —2 and Zi = 2, the redox reaction is 

Oo^^02 + 2e-+VS+. (21) 

In these reactions Vad represents a vacant oxygen ion adsorption site on top of the film and Oq 
represents an occupied oxygen site in the outermost layer of the film. We include these sites in the 
equilibrium so that the concentration of ions saturates when all sites in the relevant surface layer 
are filled. The concentrations of surface ions, 9i = [Ion], are defined so that their saturation levels 
are 9i = 1 and the concentrations of the surface sites are [lonSite] = 1 — 9i. 

One can write mass-action equilibria for these redox reactions, taking into account the external 
voltage difference between the bottom electrode and the surface (since the electrons are assumed 
to reside at the bottom electrode, while the ions reside at the surface). These are given by 

Oi r-AG°-z,.eV, 



Po. ( \t ) ' ^^^^ 

where AG° is the standard free energy of formation of the surface ion at = 1 bar and Vex = 0, 
and e is the magnitude of the electron charge. This expression is analogous to the Langmuir 
adsorption isotherm used in interfacial electrochemistr})^ for adsorption of neutral species onto a 
conducting electrode exposed to ions in a solution. Here, we consider adsorption of ions onto a 
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polar surface exposed to neutral species in a chemical environment. Thus our Vex is the potential 
of the adsorbed ions relative to the electrons, rather than the potential of the electrons relative to 
the ions in solution as in the typical electrochemical case. 

The standard free energies can depend not only on temperature but also on the polarization of 
the film, since the surface structure changes with P. For simplicity we assume that they can be all 
described by the same parameter A' using 

AG° = ACriT) + {z,eX'/eo)P. (23) 

If A' is positive, then a more positively polarized film tends to stabilize negative surface ions, and 
vice versa. Note that the effect of A' is in addition to the electrostatic energy already included 



through the ZicVex term in Eq. (22). The density of free charge on the surface is the sum of those 



from the various surface ions, giving 



(24) 



A. 



where the A- ^ are the saturation densities of the surface ions. 

Using Eqs. (22)-(24) and Eq. ([T]), the surface ion concentrations 9i can be calculated for given 



Vex and po2- Parameter values for a system with one positive surface ion i = ® and one negative 
surface ion i = Q are given in Table |II| here AC^ and AGq are taken to be independent of 
temperature. We have assumed that the saturation densities of the surface ions A^^ are both one 
per PbTiOs unit cell area. For divalent surface ions, this saturation density would provide more 
than twice the charge density needed to fully compensate the typical polarization of PbTiOa. 
Figure [TT| shows 9i as a function of Vex for several different values of po2- Changing shifts 



the external voltage scale by [kT/ {zinie)]A\npo2 for each ion. Figure 12 shows the corresponding 
surface charge density a. Three plateaus occur - two at extreme values of Vex, where one or the 
other of the ionic surface concentrations 9i saturates at unity, and a third near zero surface charge 
over the range of Vex for which both 9i are small compared to unity. From the shape of the charge 



vs. voltage curves in Fig. 12, one can see that this fixed po2 boundary condition has regions that 
correspond to fixed a separated by regions that correspond approximately to fixed Vex- Thus fixed 
P02 does not correspond to either fixed a or fixed Vex- As we shall see, this strongly affects the 
equilibrium phase diagram. 



The values of P02 and Vex that give a = can be obtained by solving Eqs. (22)-(24) for 
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a = P = 0, 6i « 1, to give 



In 



P02 



Uq - riQ 



kT 



(25) 



where is the temperature-dependent oxygen partial pressure that gives a = P = at V^x = 0, 



—nQ-n 



Uq-Uq 



^77^ ^ + ln 

kT 



(26) 



As shown below, this value of marks the transition between oppositely polarized films on the 
phase diagram. 



A. Equilibrium solutions at controlled oxygen partial pressure 

Equilibrium solutions can be calculated by obtaining a relationship between Em and P due 
to the fixed P02 chemical boundary condition, and solving it simultaneously with the constitutive 
relation for the ferroelectric, as we did for the electronic case in Fig. [2j The solution for a as a 
function of V^x shown in Fig. 12 along with Eqs. ([!]) and (|2]) gives a relation between and P 



for a given P02, which can be solved simultaneously with Eq. (|3]) to obtain the overall equilibrium. 



This is illustrated in Fig. 13 where the chemical boundary condition for three values of is 



shown. Changing shifts the boundary condition curve along the Ein axis. The boundary 



condition curve is centered on P = 0, Ein = when P02 is equal to Pq° of Eq. (25 1. 

An approximate expression for the chemical boundary condition can be obtained by making 
some simplifying assumptions. The dependence of a on V^x obtained from the depolarizing field, 
Eq. ([T]), is negligible compared with that obtained from the chemical equilibria, Eq. ( [24] ), shown 
in Fig. [12] One can approximate Eq. Q as a ^ —P. In addition, we can neglect one of the ion 
concentrations 9^ or 9q relative to the other, depending on the sign of the film polarization. This 
leads to the limiting expressions 

V;^ ^ _^ 1„ -^P'J^ _ ^ _ + ^li^, (27) 

ZiC l + AiP/ziC ZiC eo ZiUiC 

for i = or © (positive or negative film polarization, respectively). For the parameters used in 



Fig. [12] e.g. far below Tc, the approximation (27) is very close to the exact solution obtained 



numerically. Substituting the approximation (27) into Eq. ([2]), one obtains relationships between 
field and polarization given by 



ZiC l~\-AiFIZie ZiC ZiUiC ' 



eoEm > (28) 
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for 2 = or ©. Here we have introduced a new parameter, defined by 

= A - A72. 

If A' is positive then A^^ is smaller than A, and in particular A^^ can be negative. 



(29) 



B. Stability of equilibrium solutions 



Because of the plateau in the P{Ein) shape of the chemical boundary condition shown in 



Fig. 13 there can be as many as five equilibrium solutions given by the intersections. A total 



free energy function that is minimized at equilibrium can be used to determine which solutions 
are stable, metastable, and unstable. The Gibbs free energy consistent with the above treatment of 
ionic surface compensation is 

[a + Pf 



j=e,9 



kT 



2eo 
kT 



Ata^ 

+ 

rii 



(30) 



Minimizing this Q with respect to P at constant P02, 9q, and 6'© (and therefore constant a) gives 



1 dg_ 

t dP 



a + P 
eo 



(31) 



This agrees with the constitutive relation, Eq. ([3]), like the case for electronic compensation, 



Eq. (14). As in that case, because of the large value of C /Tc for PbTiOs, the equilibrium po- 
larization is given to a good approximation by P ^ — a. The free energy expression then becomes 

,2 



G ^ tf{a) + 



av 

eo 



i=e,i 



kT 



9,AGr 9,lnpo2 



kT 



Hi 



+ 9,\n9, + {l-9i)Hl-9i) 



(32) 



When the derivatives of this free energy with respect to 9^ and 9q at fixed P02 are set to zero, this 
yields the equilibrium relations given above in Eqs. (|22 24 1, where Vex of Eq. ([T]) is now given by 

Vex^2Xa/eo + tf'ia). 
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The global minimum of the free energy Q of Eq. (32 1 with respect to 9q and 9^ typically occurs 
either at ^ or 9^ ^ 0. The generality of this result can be evaluated by re-expressing the 9q 



and terms in the free energy Eq. (32 1 using new variables a and S = Zq 
Minimizing Q with respect to 5 at fixed a gives 







2e 



dg 



+ 
+ 



85 
kT 



In 



9n 



ZqUq 



kT 



In 



9, 



e 



ZmTiff^Zc^n, 



e"'e 



kT lnpo2- 



(33) 



The first term is positive in cases such as the one we consider, when there is a region of intermediate 
P02 with low concentrations of both positive and negative surface ions. The third term is zero for 
typical values of the Zi and rij. Thus the equilibrium condition requires that the second term 
be negative, which occurs only when either 9q or 9^, is very small. Substituting this result into 



Eq. (32) gives 



g ^ tf{a) 



eo 



+ 



kT 



m 



ZiC 



ZiC V 



kT lnpo2 



(34) 



ZiC 



In 1 - 



Act 



ZiC 



where i = for 6^® ^ 0, positive P, and negative cr, or i = © for 9q ~ 0, negative P, and positive 
a. The third term is proportional to cr, like a field term, but the constant of proportionality changes 
when cr changes sign and the ionic species at the surface change between positive and negative 
ions. This change in slope of g{a) at a = can produce a stable or metastable minimum. 

and 
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At intermediate 



The free energy of Eq. (34) is plotted versus P02IP02 ^'^^ ^ Figs, 
values of P02, there are three (meta-)stable solutions corresponding to local minima in g{<j), at 
positive, zero, and negative polarization. These equilibrium solutions satisfy the equations of state 







da eo ZiC \ 



kT In P02 



rii 



+ 



kT 



ZiC 



In 



Aor 



ZiC 



In 1 



A(T 



ZiC 



(35) 



and the limits of metastability of these solutions can be obtained from 

n ^'^~.f'V ^^2At kT 



Zjca 1 1 



Zi.e 



(36) 
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where i = or © as in Eq. p4\ . Figure [T6[a,b) shows the polarizations and energies of these 
solutions as a function of poi/Poi - energy of the solution at P = is zero, while the energies 
depend on Po2lP°02 other two solutions. For the parameters used here, e.g. a 3.2 nm film 

thickness, the energies of all three solutions are almost equal at P02IP02 ~ ^- "^^^ P = solution 
will be the stable (global minimum) solution for thinner films at intermediate po2- Here this 
solution is stable against a nonuniformity, unlike the electronic compensation case. For example, 



Fig. 16 c,d) shows the results for a 1 .6 nm thick film, with all other parameters the same. In such 
thin films, where the central fiat region of the boundary condition P(£'j„) curve spans a large 
range of i?j„, the positive and negative solutions do not overlap, and the P = solution is the only 
solution for the range of P02 where both 9^ and 9q are small. 

For the polar phases, the last term of Eq. ([36]) is typically small enough, except near Tc, that this 



condition for the instability is very similar to those for the electronic compensation cases, Eqs. ( 12 ) 



and (18). Thus at the metastability limit of the polar phases, the internal field reaches the same 
intrinsic coercive field in the ionic compensation case as it does in the electronic compensation 
cases. 



C. Phase diagram for controlled P02 

The effect of ionic surface compensation on the ferroelectric phase transition can be explored 
by solving for the polarization and field as a function of temperature as well as P02 and film 
thickness. As can be guessed from the fixed-temperature results shown above, the temperature 
dependences of P, Em, and the Curie point Tc (i.e. the temperature of the equilibrium boundary 
between the polar and nonpolar phases) all vary with the P02 of the environment. 



Figure 17 shows equilibrium polarization phase diagrams as a function of T and 
various film thicknesses. In addition to the stable and metastable equilibrium phase boundaries, 
the metastability limits of the polar and nonpolar phases are shown. These phase diagrams are cal- 
culated using parameter values given in Table |llj The oxygen pressure scale has been normalized 
topo° (T), which produces symmetric diagrams when Uq = — n®, Zq = —z^. 



The equilibrium phase diagrams as a function of P02 for ionic compensation. Fig. 17, differ 
qualitatively from the standard second-order ferroelectric phase diagrams as a function of Vex or 
a for electronic compensation, Figs.[5]and|9j The ionic phase diagrams show temperature ranges 
where the nonpolar phase is stable at intermediate P02 separating the positive and negative polar 
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phases at high and low poa, respectively. As film thickness becomes smaller, this "wedge" of 
nonpolar phase extends to lower temperature, reaching K for thicknesses less than about 1 nm 
for the parameter values used here. For films with smaller thickness, an inverted ferroelectric 
transition remains at extreme values of poj, with the polar phase stable at temperatures above the 
phase boundary, and the nonpolar phase stable below the boundary. For thicker films, there is a 
triple point where the first order transitions between the positive and negative polar and nonpolar 
phases meet at P02/P02 ~ while at extreme values of P02 there is no phase transition as a 
function of T between the polar and nonpolar phases, similar to the case at nonzero Vex in Fig.|5]for 
electronic compensation. For all film thicknesses, the high temperature ends of the polar/nonpolar 
phase boundaries terminate at two critical points. In a range of temperature below these critical 
points, the regions of (meta) stability of the positive and negative polar phases do not overlap. Here, 
switching transitions between oppositely polarized states at fixed T driven by changing must 
occur through an intermediate nonpolar state. 

The appearance of the nonpolar phase between the polar phases at lower temperature is directly 
related to the non-linear dependence of surface charge a on ln(j9o2) fixed V^x- This has a plateau 
at a value near cr = for intermediate po^ values, where the concentrations of both positive and 
negative surface ions are small. The low value of a in this region can be insufficient to stabilize 
either polar phase. 



Figure 18 shows the temperatures as a function of film thickness of the critical points, the 
polar phase instabilities at ~ triple point. The triple point is the minimum 

equilibrium Tc- At temperatures between the triple point and the critical points, the nonpolar 
phase intervenes between the polar phases at equilibrium. An expression for the temperatures Tcr 
of the critical points can be obtained by setting the second and third derivatives of the free energy 
simultaneously to zero, Eq. (36) and 

d^g kT{i-^-jf 

Zita^ f 1 - ^ 



° = a^ — ^^^^ 



where i = or © as in Eq. ( [34| ). In the approximation that a is small at the critical point, these 
reduce to 

0..(a5(T„,.^)%2„.3(|^)\ (38) 
We use the double ~~ symbol to indicate a rough approximation, in this case because it becomes 



invalid at small t. Nonetheless, Eq. p8| ) shows that the temperatures of the critical points for the 
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ionic compensation case are suppressed by an additional thickness-dependent term not present in 
the electronic compensation case, Eq. ( 13 1. Even if the effective screening length is zero, = 0, 
the Tcr are changed by an amount 

AT,, ^ - T° -2eoC{2al,y/' ( 4^ ) ' ^^^^ 
Using the LGD parameterPfor PbTiOa coherently strained to SrTiOs, = 2, and a thickness of 
t = 3.2 nm, one obtains ATcr ~~ — 115 K. 



Figure 19 shows the internal field (along with the phase boundaries) as a function of P02/P02 



and T for a 1.6 nm thick film. Like the electronic compensation cases. Figs. [6] and 10 the in 



ternal electric field is inverted in the polar phases near the phase boundaries because of the in- 
completely neutralized depolarizing field. While the electronic compensation model requires a 
nonzero screening length A to produce an inverted field, the ionic compensation model does not. 
The magnitude of the inverted field at the phase boundary is much larger for ionic than for elec- 
tronic compensation in the cases shown. In all cases the inverted field regions extend above the 
critical point(s). The oxygen partial pressure corresponding to zero internal field can be obtained 
by setting the numerator in Eq. ( [28] ) to zero, giving 

1 E,„=o _ 2niZieX^Po f -APo/zie \ 
m Pn ^ — h rij In — 

where Pq is the T-dependent zero-field spontaneous polarization of the epitaxially strained film 
given by the solution to Eq. ([2]) with Ein = 0, and i = or © for positive or negative values of 
Pq. As for electronic compensation, the conditions for zero field are independent of film thickness, 
and they intersect at T^. Rough values of the oxygen partial pressure at the critical points Pq^ can 



be obtained by assuming that the field is zero and neglecting the first two terms in Eq. (40). This 
gives kTcr Inp^^ ~~ ?T,jAG°°, where z = for positive P (high i = for negative P 

(low P02). 

D. Phase diagram for controlled surface oxygen density 

It is instructive to plot equilibrium phase diagrams for ionic compensation as a function of the 
net excess surface oxygen density, defined by 
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In the typical case where we can neglect one or the other of the 9i for positive or negative a, one 
obtains po ~ 2a /riiZie, for z = or ©, so that po is simply proportional to a for each range. 



Figures 20 and 21 show the polarization and internal field plotted as a color scale on po vs. T axes 
for t = 1.6 nm. These correspond to the vs. T diagrams shown in Figs. [T7|^b) and 19 Such 
controlled po diagrams can be obtained using the Helmholtz free energy 



A = g + P^kT\npo, 



eo 



(42) 



+ 



kT 



In 



1 - 



A, 



a 



ZiC 



In 1 



A, 



a 



ZiC 



where z = or © as in Eq. 34 



Because po is a conserved order parameter, the phase boundaries occuring on the controlled 



potential diagrams become two-phase fields on Figs. 20 and 21 as in the electronic compensation 



case for controlled surface charge, Figs. |9] and [TOj These show the phase separation that would 
occur in a closed system if the net amount of excess surface oxygen is fixed, rather than the 
external P02 ■ The phase diagram for ionic compensation is more complex than that for electronic 
compensation; there are two equilibrium two-phase fields between polar and nonpolar phases 
above a tie-line at the triple point temperature. 



V. DISCUSSION 



The new parameters in the model developed above for ionic surface compensation are n^, Zi, Ai, 
and AG°° for i = Q and 0, as well as A'. These can be related to the locations of the features on 
the phase diagram. Approximate expressions are given above that show how the AG°° determine 
and the Pq^, which give the center and width of the phase diagram in P02 coordinates. These 
expressions are particularly simple when the phase diagram is symmetric in coordinates of P02/P02 
versus T, i.e. when = —Uq, Zq = —z^, and Aq = A^. In this case one obtains Inp^" = 
ne(AGg' - AQ°)/(2A;T), ln(p- ) ±nQ{AG°^ + AG°°)/{2kT,r). The value of A' 
affects the critical temperatures Tcr, which may be suppressed or enhanced relative to the electronic 
compensation value. 

Although the phase diagrams we have shown are symmetric when plotted in P02/P02 versus T 
coordinates, the value of Pq^ is expected to be a function of T. Thus experimental phase diagrams 
obtained as a function of P02 versus T are not expected to be symmetric. Figure [22] shows the 
trajectories of constant P02 on the same T vs. P02/P02 ^^^^ "^^^ '^^ P^*^*- phase diagrams. 
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Note that we have neglected any temperature dependence of AG°° in these calculations. Such 
temperature dependence might be expected since the entropy of O2 in the environment may be 
different than that of the adsorbed ions. 

While the Gibbs free energy expressions for the ionic and electronic compensation cases are 
very similar, the phase diagrams differ qualitatively because a different parameter is fixed. If we 
neglect the polarization dependence of the AG° so that A' = and we consider 9i « 1 for both 



i = Q and ©, then by using the mass action equilibria Eqs. (22 1 one can show that the Gibbs free 
energy for ionic compensation, Eq. (30), reduces to Eq. (|9]) used for the fixed Vex case. However, 
the fixed Vex and fixed P02 conditions lead to different equilibrium free energy surfaces. Figs. |3] 



and 14 even in the case of A' = and 6i « 1. The relationship between the fixed Vex and 



fixed conditions can be seen from Eq. (27 1. Here enters into the expression for Vex simply 
through the term kTlnpo^/iziUie). In general, however, fixed does not correspond to fixed 
Vex because there are other terms and they depend upon polarization. In particular, if the value of 
AG°°/ (zjc) differs for positive and negative surface ions, then there is an abrupt jump in Vex when 
crossing P = 0. Thus fixed Vex and fixed po^ constraints produce different equilibrium behavior 
even when the free energy expression for both cases is the same. 



The form of the phase diagrams in Fig. 17 in which a stable nonpolar phase intervenes between 
the polar phases is due to the appearance of a third local minimum in Q (cr) near cr = as shown in 



Fig. 15 This is conceptually similar to the behavior of a ferroelectric with a first-order transition, 
for which the coefficient of in the free energy expression is negative.'^^'^ For example. Figs. 23 



and |24] in Appendix |A] show the phase diagrams for controlled internal field Ein and net surface 
charge a for unstressed bulk PbTiOs with ideal electrodes (A = 0). Here the topology of the 



equilibrium phase boundaries is similar to that in Figs.[T7[a,b) and [20} with a triple point and two 
critical points. However, the range of temperatures spanned by this structure in ultrathin films with 
ionic compensation can be much larger than in bulk PbTiOa. Furthermore, the polarization of the 
paraelectric phase at temperatures below the critical points is much closer to zero in ultrathin films 
with ionic compensation, because of the sharp minimum in ^(cr) at a = 0. 

The appearance of a stable nonpolar state between the polar states on the P02 vs. T phase 
diagram can affect the mechanism of switching and the internal field at which switching occurs 
(i.e. the coercive field). During switching by ramping P02, the film may first become unstable with 
respect to the nonpolar state before reaching P02 values that stabilize the opposite polarization, thus 
suppressing nucleation of oppositely polarized domains. In this case the internal field could reach 
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the intrinsic coercive field and switching occur by a continuous, spinodal mechanism without 
nucleation. This could produce the recently observed crossover to a continuous mechanisrcP^ 
through an equilibrium pathway not requiring kinetic suppression of nucleation. 

The model developed here contains several assumptions that could be relaxed in future exten- 
sions. We assume that the effective screening length A is not negative, so that electronic interfacial 
effects tend to suppress rather than enhance polarization in ultrathin films. We also neglect any 
polarization dependence of A. Ab initio calculation^^'^ indicate that in some systems the inter- 
faces enhance film polarization, which can be modeled with a negative A, and that A depends on 
P. These effects could be included by modifications to our electrostatic boundary conditions and 
free energy expressions. We constrain the free and bound charge at each interface to reside in 
single planes, so that there is no space charge. Such space charge could be included as has been 
done previously in models with semiconducting ferroelectric films and/or electrodes .'^'^'^ Since 
the screening layer of thickness A is a conceptual construct rather than an actual dielectric layer 
in our model, we do not consider tunneling of free charge across this layer, which has been re- 
cently considered for systems with a dielectric separating the electrode from the ferroelectric.'^ 
These effects could be added for such systems. We also neglect the possibility of equilibrium 180° 
stripe domain formatiorP^ in which nanoscale domain structures reduce the depolarizing field 
even when there is little or no electronic or ionic compensation charge at one or both interfaces. A 
full treatment of equilibrium stripe domains for the ionic compensation case would be valuable in 
future work. 



VI. SUMMARY AND CONCLUSIONS 



Ionic compensation of a ferroelectric surface due to chemical equilibrium with an environment 



introduces new features into the phase diagrams. Figs. 17 andEOl not present in the standard phase 



diagrams for a second-order transition in a film with electronic compensation. Figs. |5] and [9j The 



constant P02 chemical boundary condition shown in Fig. 13 is a hybrid between the constant Vex 
and constant a boundary conditions shown in Fig. |2j Because the surface concentrations of ionic 
species 9i are limited to values between zero and unity, constant surface charge regimes occur 
when the 9i are saturated. In the regimes where one of the 9i is varying between these limits, the 
boundary condition is similar to a fixed Vex condition. There are two independent relations for the 
surface charge cr as a function of P02, depending upon whether positive or negative surface ions 



22 



predominate. In the region where there is insufficient surface charge of either sign to stabilize 
a polar state, the nonpolar state becomes stable between the positive and negative polar states, 
producing two critical points, a triple point, and a strong dependence of Tc on po^- Large inverted 
internal fields occur at equilibrium in the polar phases near the phase boundaries with the nonpolar 
phase. Manipulation of ultrathin ferroelectric films via controlled ionic compensation may thus 
allow experimental access to exotic nonpolar and high-field states such as those modeled in recent 
ab initio calculation^^^^^ that would not be stable under electronic compensation conditions. 
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Appendix A: Phase Diagrams for First-Order Ferroelectric 

The phase diagrams for a bulk ferroelectric that has a first-order transition at zero field are 
similar to those for ultrathin films with ionic surface compensation. Here we present calculated 
phase diagrams showing the region near the critical points for the weakly first-order transition in 
unstressed bulk PbTiOa. The features in these phase diagrams can be compared with those for ul- 
trathin epitaxially strained films shown above, which have a second-order transition for electronic 
compensation but a strongly first-order transition for ionic surface compensation. 

These phase diagrams are calculated using the Landau-Ginzburg-Devonshire expression for the 
free energy per unit volume of unstressed bulk PbTiOs with ideal electrodes having zero screening 
length (A = 0), 

= aiP^ + auP^ + «inP^ + eoEl/2 - E^^P (Al) 



with «! = (T — To)/2eoC, where the parameters are given in Table I of the main paper. Figure 23 



shows the phase diagram for controlled internal field and Fig. 24 shows the phase diagram for 
controlled net surface charge a. These are typical for a ferroelectric with a first-order transition, 
for which the coefficient of P'^ in the free energy expression is negative. The topology of the 
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equilibrium phase boundaries is similar to the ionic compensation case, Figs.[T7[a,b) and|20| with 
a triple point and two critical points. However, the paraelectric phase above the triple point has a 
much wider range of polarization around zero in this case, reflecting the relatively broad minimum 
in near P = 0. 
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FIG. 1: Schematic of polarization, displacement, electric field, and electric potential in the bulk and at the 
interfaces of a ferroelectric film of thickness t and polarization P. Compensating planes of charge density 
a can be considered to reside at a separation A equal to the effective screening length in the electrodes. 
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FIG. 2: (Color online) Polarization vs. internal field. Straight red lines show (a) fixed Vex (at t = 3.2 nm 
and X/t = 10^^) or (b) fixed a boundary conditions from Eq. (j2]l, for the three values of the fixed quantity 
given in the legend. In each case, the "S" shaped blue curve is the constitutive relation, Eq. ([3]), for PbTiOs 
coherently strained to SrTiOs at 644 K. Marked intersections correspond to stable or metastable equilibrium 
solutions. 
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FIG. 3: (Color online) Gibbs free energy Q vs. P and Vex for a t = 3.2 nm PbTiOs film coherently strained 
to SrTiOs at 644 K, with X/t = 10"^. Color scale gives values of Q. 
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FIG. 4: (Color online) Equilibrium solutions for (a) polarization P, (b) strain ^3, and (c) Gibbs free energy 
^ as a function of Vex, calculated for PbTiOs coherently strained to SrTiOs at 644 K with t = 3.2 nm 
and X/t = 10"^. Red (blue) curves are positive (negative) polarization; solid (dashed) segments are stable 
(metastable); closed (open) circles show equilibrium transition (instability) points. 
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FIG. 5: (Color online) Equilibrium polarization phase diagram as a function of Vex and T for PbTiOs 
coherently strained to SrTiOs with t = 3.2 nm and X/t = 10^^. Color scale gives polarization of stable 
phase. Solid black line is phase boundary between positive and negative polar ferroelectric (F-i- and F-) 
phases, terminating at Tq. Dashed red and blue curves are metastability limits of the F-i- and F- phases, 
respectively. 
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FIG. 6: (Color online) Internal field in region near the Curie point, corresponding to Fig. [5] Color scale 
gives electric field in stable phase. Dashed black curves show conditions for zero field, which intersect at 
(open circle), the Curie point for a film without depolarizing field (A = 0). 
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FIG. 7: (Color online) Helmholtz free energy A vs. polarization P and surface charge o" for a t = 3.2 nm 
PbTiOa film coherently strained to SrTiOs at 644 K, with X/t = 10~^. Color scale gives values of A. To 
emphasize the equilibrium solutions, only the region near P = —a is plotted since A is very large outside 
this region. 
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FIG. 8: (Color online) Equilibrium solutions for (a) external voltage Vex and (b) Helmholtz free energy 
^ as a function of surface charge a corresponding to Fig. [7] Red (blue) curves are positive (negative) 
polarization; all values are stable with respect to P variations when a is spatially uniform. Closed (open) 
circles show the equilibrium transition (instability) points when a can be nonuniform. 
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FIG. 9: (Color online) Equilibrium phase diagram as a function of net surface charge a and T when a can 
be nonuniform, for PbTiOs coherently strained to SrTiOs with t = 3.2 nm and \/t = 10~^. Color scale 
gives polarization in single-phase region. Solid black line is phase boundary between positive and negative 
polar ferroelectric (F+ and F-) phases and a two-phase field, which terminates at Tc (filled circle). Dashed 
red and blue curves are metastability limits of the F-i- and F- phases, respectively. 
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FIG. 10: (Color online) Internal field in the region near the Curie point, corresponding to Fig.[9] Color scale 
gives electric field in single-phase region. Dashed black curves show conditions for zero field. 



36 



10" 



10 



o 
U 



10 







\ s V / / 
\ \/^ f f 


lo--^ 

— 10-6 

---10-^ 


// \\\ 

/ \ V 



0.6 0.4 0.2 -0.2 -0.4 
External Voltage (V) 



FIG. 11: (Color online) Concentrations of positive 00 (red curves) and negative Qq (blue curves) surface 
ions as a function of calculated for the values of -po^ (bar) given in the legend at T = 644 K, t = 3.2 nm. 
Parameter values used are given in Table [ill 
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FIG. 12: Surface charge density cr as a function of Vex at three values of po2, corresponding to the concen- 
trations shown in Fig.[TT] Parameter values used are given in Table 
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FIG. 13: (Color online) Polarization vs. internal field relationships arising from the constitutive relation 
for the ferroelectric film (blue "S" shaped curve) and from the chemical boundary condition (red curves 
with plateau) for values shown in legend (bar). Marked intersections correspond to stable or metastable 
equilibrium solutions. Parameter values used are given in Tables [I] and [IH with T = 644 K, t = 3.2 nm. 
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FIG. 14: (Color online) Free energy as a function of a and VOilVo^- P^^mcter values used are given in 
Tables [l] and [nj with T = 644 K, t = 3.2 nm. Color scale gives values of Q. 
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FIG. 15: (Color online) Free energy as a function of a at three given values of po2/Po2- P^^^i^tsi" values 
used are given in Tables [l| and [lH with T = 644 K, t = 3.2 nm. 
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FIG. 16: (Color online) Polarization P and free energy Q of the (meta-)stable equilibrium solutions as a 
function of P02/Po2 - ^1"^' green, and red curves are for negative, zero, and positive polarization solutions; 
solid and dashed regions are stable and metastable, respectively. Parameter values used are given in Tables |l] 
and|Tl| with T = 644 K, and t = 3.2 nm for plots (a) and (b), t = 1.6 nm for plots (c) and (d). 
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FIG. 17: (Color online) Phase diagrams as a function of "Po^lVo^ ^'^'^ ^ f^J" three thicknesses of PbTiOs 
coherently strained to SrTiOs, using parameters in Tables [l] and [Tl] Color scale gives polarization. Solid and 
dash-dot black curves are stable and metastable phase boundaries, respectively, between nonpolar paraelec- 
tric (?) and positive and negative polar ferroelectric (F+ and F-) phases. Dashed red, blue, and white curves 
are metastability limits of the F+, F-, and P phases, respectively. 
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FIG. 18: The temperatures of the critical points T^r, the polar phase instabilities at = Pq^, and the 
triple point, as a function of film thickness. 
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FIG. 19: (Color online) Internal field for t = 1.6 nm, corresponding to Fig.[T7jb). Color scale gives electric 
field in stable phase. Dashed black curves show conditions for zero field, which intersect at T^. 
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FIG. 20: (Color online) Phase diagram as a function of net excess surface oxygen density po and T for 1.6 
nm of PbTiOa coherently strained to SrTiOa, corresponding to Fig. [ITjb). Color scale gives polarization 
in single-phase region; three two-phase regions are marked. Solid and dash-dot black curves are stable 
and metastable phase boundaries, respectively, between nonpolar paraelectric (?) and positive and negative 
polar ferroelectric (F-i- and F-) phases. Dashed red, blue, and black curves are metastability limits of the F-i-, 
F-, and P phases, respectively. 
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FIG. 21: (Color online) Internal field corresponding to Fig. 20 Color scale gives electric field in stable 
phase. Dashed black curves show conditions for zero field, which intersect at (open circle). 
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FIG. 22: Curves show fixed values of (bar) given in the legend, plotted using same normalized axes 
used for the phase diagrams, with parameters corresponding to Table |ll] 
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FIG. 23: (Color online) Region near critical points on equilibrium polarization phase diagram as a function 
of Ein and T, for unstressed bulk PbTiOs with ideal electrodes (A = 0). Color scale gives polarization 
of stable phase. Solid and dash-dot black curves are stable and metastable phase boundaries, respectively, 
between nonpolar paraelectric (P) and positive and negative polar ferroelectric (F-i- and F-) phases. Dashed 
red, blue, and white curves are metastability limits of the F-i-, F-, and P phases, respectively. 
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FIG. 24: (Color online) Region near critical points on equilibrium polarization phase diagram for unstressed 
bulk PbTiOa with ideal electrodes (A = 0) corresponding to Fig.|23} here plotted as a function of net surface 
charge a and T. Color scale gives polarization in single-phase region, which is separated by phase boundary 
(solid black curve) from three two-phase fields; metastable phase boundaries are dash-dot curves, and two 
critical points are marked by filled circles. Dashed red, blue, and black curves are metastability limits of the 
F-H, F-, and P phases, respectively. 
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Tables 
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TABLE I: Material parameterP^H] for PbTiOs. 

To 752.0 (K) Qii 8.9 x lO^^ (m'^/C^) 

C 1.5 X 10^ (K) Qi2 -2.6 X lO^^ (m'^/C^) 

an -7.25 x lO"^ (Vm^/C^) sn 8.0 x W~^^ (m^/N) 

am 2.61 x 10^ (Vm^/C^) si2 -2.5 x lO^^^ (m^/N) 
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TABLE n: Values of ionic surface compensation coefficients used in displayed plots. 



\G°^ 1.00 (eV) AG°^ 0.00 (eV) 

—2 n© 2 

2 zq —2 

A(s 1.6 X 10-19 (m2) Aq 1.6 X IQ-^^ (m^) 

A (m) A' (m) 
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